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Abstract: 

We discuss a general framework for the inclusion of heavy quark mass contributions to 
deep-inelastic structure functions and their perturbative matching to structure functions 
computed in variable-mass schemes. Our approach is based on the so-called FONLL 
method, previously introduced and applied to heavy quark hadroproduction and pho- 
toproduction. We define our framework, provide expressions up to second order in the 
strong coupling, and use them to construct matched expressions for structure functions 
up to NNLO. After checking explicitly the consistency of our results, we perform a study 
of the phenomenological impact of heavy quark terms, and compare results obtained at 
various perturbative orders, and with various prescriptions for the treatment of subleading 
terms, specifically those related to threshold behaviour. We also consider the heavy quark 
structure function F^ c and discuss issues related to the presence of mass singularities in 
their coefficient functions. 
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1 Introduction 



Interest in the inclusion of heavy flavour contributions to deep-inelastic electroproduction 
structure functions was recently revived by the discovery pQ that mass-suppressed terms 
in global parton fits can affect predictions for the total W and Z production at the LHC 
by almost 10 %. Indeed, until very recently, many global studies of parton distributions 
(PDFs, henceforth) were performed assuming heavy quarks to decouple from the structure 
functions for Q 2 < m?, but to be massless when Q 2 > m^, despite the availability of more 
refined approaches where contributions with full mass dependence are consistently taken 
into account. 

A technique for the inclusion of heavy mass-suppressed contributions to structure 
functions was developed long ago [21[3], based upon a renormalization scheme with explicit 
heavy quark decoupling [3J. Several variants of this method (usually called ACOT) were 
subsequently proposed, such as S-ACOT [5] and ACOT-x [SlE]- However, the ACOT 
method was first used for an actual global parton fit only recently, in Refs. [HE]. An 
alternative method (sometimes called TR) has also been advocated [5HT0]. and used for 
parton fits [TTHT3] . Recently, however, the methods used by the CTEQ p] and MSTW [HJ 
groups for their current parton fits, based respectively on the ACOT [2j,[3] and TR [HHTU] 
procedures, have adopted at least in part a common framework: they have been compared 
recently in Refs. [HICE], thereby elucidating differences and common aspects. 

A somewhat different technique for the inclusion of heavy quark effects, the so called 
FONLL method, was introduced in Ref. [T7] in the context of hadroproduction of heavy 
quarks. The FONLL method only relies on standard QCD factorization and calculations 
with massive quarks in the decoupling scheme of Ref. [4j and with massless quarks in the 
MS scheme. The name FONLL is motivated by the fact that the method was originally 
used to combine a fixed (second) order calculation with a next-to-leading log one; however, 
the method is entirely general, and it can be used to combine consistently a fixed order 
with a resummed calculation to any order of either. 

It is the purpose of this paper to present the application of the FONLL scheme to 
deep-inelastic structure functions. We shall see that, thanks to its simplicity, the method 
actually provides a framework for understanding differences between other existing ap- 
proaches, and for a study of the effect of different choices in the inclusion of subleading 
terms. We shall use this scheme to perform detailed studies of the phenomenological im- 
pact of the treatment of heavy quarks up to NNLO, and of the ambiguities involved in 
the procedure. 

In section 2 we will describe in detail the theory of the FONLL method for structure 
functions. In section 3 we will provide explicit expressions up to NNLO, and discuss issues 
related to perturbative ordering. In section 4 we shall compare the structure functions 
determined with different heavy quark matching procedures at various orders, and with 
various treatments of threshold subleading terms. Finally, in section 5 we shall discuss 
issues related to the presence of mass singularities in the definition of the heavy quark 
structure functions i^c an d i^b- F° r ease of reference and for definiteness, the complete 
explicit FONLL expressions for the structure function F2 are collected in an Appendix. 
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2 The FONLL method in deep- inelastic scattering 



The basic problem in the treatment of heavy quarks in a QCD process stems from the 
fact that QCD calculations are usually performed in a decoupling scheme [1] , rather than 
in the MS scheme. Indeed, in the standard MS scheme, heavy quarks contributions are 
present at all scales: for instance, the /3 function depends on nj with nj = 6 at all scales. 
In a decoupling scheme, instead, in the computation of a process characterized by the 
hard scale Q 2 all quarks with mass m 2 > Q 2 ("heavy", henceforth) are not treated in 
the MS scheme. Rather, 'heavy flavour graphs' are subtracted at zero momentum0 The 
important consequence of this definition is that heavy quarks decouple for scales much 
lower than the heavy quark mass. This implies that the GLAP evolution equations, and 
the running of a s are identical to those which one would get in the MS scheme, but with 
rij equal to the number of light flavours, nt = ni . If for Q 2 < m| one neglects all 
terms suppressed by powers of Q 2 /m^, one obtains a so-called zero-mass scheme, where 
all quarks are treated as massless, but heavy quarks are absent at Q 2 < m\. This scheme 
can be combined with the usual MS scheme for ni + 1 flavours, including the heavy quark, 
when Q 2 > m\\ if terms suppressed by powers of m\jQ 2 are neglected throughout this 
yields the so called zero-mass variable-flavour number scheme (ZMVFN from now on) of 
Rcf. [3]1 

m 2 

The ZMVFN scheme is not accurate near the threshold region, where ~ 1. The 
problem in this region is easily remedied by simply using the decoupling scheme with 
rij = nj, but retaining explicitly the full dependence on the heavy quark mass in the 
computation of hard cross sections. This way of computing however looses accuracy in 
comparison to the previous ZMVFN when Q 2 3> so that L = lnQ 2 /m 2 i 3> 1, because, 
in the decoupling scheme with rif = ni, these large logs are only included to fixed order 
in a s while in the ZMVFN scheme they are resummed to all orders. 

Heavy quark schemes are all based on the idea of matching these two different ways of 
calculating, each of which is more accurate in some kinematic region. The basic idea in the 
ACOT [3] scheme is to retain explicitly the mass dependence in Wilson coefficients of the 
ZMVFN scheme calculation (based on the massive quark factorization theorem [2] ) , while 
the basic idea in the TR |8l[9] scheme is to require continuity of physical observables in 
the threshold region, where one switches between the decoupling scheme and the massless 
calculation with nj = n\ + 1 flavours. In both cases, the matching conditions ensue from 
the requirement that computations of the same observable within different renormalization 
schemes give the same answer within the respective accuracy. 

The FONLL scheme is instead simply based on the idea of combining the decoupling 
scheme computation with the ZMVFN computation, and subtracting double counting 
terms between the two order by order in an expansion in powers of a s and L. Effectively, 

1 Here 'heavy flavour graphs' are defined [4] (recursively) as either graphs which contain a heavy flavour 
line, or counterterms to heavy flavour graphs. 

2 Note that here and henceforth we call "variable-flavour number" scheme (VFN, and thus also ZMVFN) 
a scheme where all large Q 2 logs are resummed, so in particular, when Q 2 >> m| all In terms are 
resummed. In contrast, in some previous references, VFN or ZMVFN is used to denote schemes in which 
In are not resummed to all orders: in particular in Ref. 1181 VFN denotes the scheme in which logs 
are not resummed, while that in which logs are resummed is called "PDF". In Ref. [19] ZMVFN is used 
despite the fact that In are never resummed to all orders in this reference |20] , 
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this means that for Q 2 > m| one performs the calculation in the massless scheme, but 
then one replaces all terms whose mass-dependence is known with their exact massive 
expression. This can then be done to any desired order. 



2.1 FONLL expressions for DIS structure functions 

Let us now see explicitly how the FONLL scheme works in the case of a generic deep- 
inelastic structure function F(x,Q 2 ): in Sects. [3][3] we shall specifically discuss both F2 
and Fl, but for the time being the distinction is irrelevant, so we will refer to a generic 
structure function F. For the rest of the paper we will assume ni light flavours, with 
a single heavy flavour of mass = m. In Sec. HI when discussing phenomenological 
implications, we will study specifically the case of charm. 

The expression of F(x, Q 2 ) in the fully massless scheme, which is accurate when W ^> 
4m 2 , with W = yjQ 2 {\ — x)/x, is factorized in terms of PDFs^l 

f^\x,q 2 ) = x f 1 ^ £ ct' +1) (-,a^HQ 2 ))ft l+1 \y,Q 2 ), (1) 
j x y . _ , T \y j 

i=q,q,n,h,g 

where q is any light quark, h is the heavy quark, and n\ + 1 is the total number of flavours. 
In this expressions, the PDFs include the heavy flavour as a light parton, and the strong 
coupling as and the PDFs satisfy standard GLAP equations with rif = n\ + 1. 
Henceforth, we shall refer to Eq. ([TJ as the determination of F(x,Q 2 ) in the massless 
scheme. 

The expression of F(x,Q 2 ) in the decoupling scheme with nj = n;, which is accurate 
when W ~ 4m 2 , is instead given by 

f^(x,q 2 ) = x f 1 * c^ l) (~X 2 M? l) {Q 2 )) ft\y,Q 2 ). (2) 

Jx y i=q,q,g \y m J 

The coefficient functions are computed fully retaining the mass dependence, while now 
atp 1 '* and f^ ni \y,Q 2 ) obey standard MS evolution equations with ni flavours. When 

W 2 <C 4m 2 the heavy quark mass dependence drops out of the coefficient functions c\ ni \ 
which then reduce to the standard massless MS coefficient functions with n\ flavours. 
Henceforth, we shall refer to Eq. ([2]) as the determination of F(x,Q 2 ) in the massive 
scheme. 

In order to carry out the FONLL procedure, we need to express the decoupling scheme 
cross section, Eq. ([2]) in terms of ai n ' +1 ^ and fi ni+1 ^ for i ^ h,h. The coupling constant 
and PDFs are related in the two schemes by equations of the form 

00 

ai^ +1 )(Q 2 ) = a^(Q 2 ) + J>(L)x (a^(m 2 ))\ (3) 

i=2 

ft l+1 \x,Q 2 ) = f x d j- £ K l3 ^LMr i) (Q 2 ))f ( r\y,Q 2 ), W 



j=q,Q,g 



Note that we define PDFs f(x) in such a way that xf(x) is a momentum distribution; a different 
definition such that f(x) is the momentum distribution itself, which differs from our own by a factor of x, 
has also been sometimes used (see e.g. Ref. [21] ). 
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where 

L = logQ 2 /™ 2 - (5) 

The coefficients Cj(L) are polynomials in L, and the functions Kij can be expressed as an 
expansion in powers of a s , with coefficients that are polynomials in L. 

The first 2ni + 1 equations in Eq. ([3]) relate the light quark and gluon PDFs in the 
two schemes, and can be inverted to express the massive-scheme PDFs in terms of the 
massless-scheme ones. The last two equations in Eq. (j4]) express heavy quark PDFs in 
terms of the light flavour ones, under the assumption that the heavy flavour PDF is 
generated perturbatively. A possible intrinsic heavy flavour contribution |22U23j could be 
added as a separate contribution on the right-hand side of Eq. ([3]) for i = h, h. 

Inverting Eqs. (|3]|3J) and substituting in Eq. ([2]), one can obtain an expression of 
F^ n '\x,Q 2 ) in terms of ai n ' +1) and f n ' +1 ) : 



F^ ni \x,Q 2 ) = x f 1 — £ B t ( x -X 2 ,^ +1) {Q 2 )) ft +1 \y,Q% 
Jx y \y ™ 2 / 



(6) 



where the coefficient functions B t are such that substituting the matching relations Eqs. ([3} 
H]) in Eq.([6]) one gets back the original expression Eq. ([2]). We can thus use for F^ 11 ^ the 
expression given in Eq. (|6|), and avoid any further reference to cvg and n ■ from now 
on, we shall use Eq. ([6]) as the determination of F in the massive scheme. 

In order to match the two expressions for F in the massless scheme, Eq. ([I]), and 
in the massive scheme, Eq. 0, we now work out their perturbative expansion. Using 
GLAP evolution in the absence of intrinsic heavy quark contributions, the heavy quark 
PDFs fl n) , at the scale Q 2 which appear in the massless-scheme expression Eq. dU) 

can be determined in terms of the light-quark PDFs, /; with i 7^ h, h at the scale m, 
convoluted with coefficient functions which can be expressed as a power series in as \m), 
with coefficients that are polynomials in L, or, alternatively, in terms of the light-quark 
parton distributions f^ ni+1 ^ at the scale Q 2 convoluted with coefficient functions expressed 
as a power series in as™ !+1 ^((5 2 ), again with (different) coefficients that are polynomials in 
L. 

Thus, the massless-scheme expression Eq. ([1]) may be written entirely in terms of 
light-quark PDFs: 



F^(x,Q 2 ) = x £ Af l+1) (-,L,at +1 \Q 2 )) ft l+1 \y,Q% 



(7) 



i=q,q,g 



where the A^ ni+1 ^ coefficient functions are given by a perturbative expansion of the form 

N 00 , 

A^ +1) (z,L,a^ +1 HQ 2 )) = £ (^(Q^Y^Afiz) (af^{Q 2 )L) , (8) 

p=0 k=0 

where at leading order N = 0, and at N fc LO N = k. 

On the other hand, the massive-scheme expression Eq. ([6]) is also written in terms of 
the light quark PDFs, with coefficient functions Bi which admit a fixed order expansion 
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of the form 

where P is the order of the expansion needed to reach the desired accuracy. It follows that 
the sum B^' p of all contributions to the massive-scheme expression Eq. ([9]) which do not 
vanish when Q 2 3> m 2 must also be present in the massless-scheme result, i.e. they must 
correspond to the logarithmic and constant terms in Eq. ([8]). Namely, to p-th order in 

a ( r i+1 \Q 2 ) 

Bf^U^)^±^\ X )L\ (10) 
^ ' k=0 

where we denote by fi. ' p the massless limit of Bf , in the sense that 



lim 

m— ¥0 



2\ , , / n2 



0. (11) 



In other words, B^ ,p is obtained from Bf by retaining all logarithmic and constant terms, 
and dropping all terms suppressed by powers of m/Q. 

The FONLL method can be simply stated as follows: in the massless-scheme expres- 
sion Eq. ([7]) one replaces all contributions to the expansion of the coefficient functions 

(z, L, a ( P l+1) (Q 2 )^ which appear in Bf ] ' p (x, , Eq. CED, with their fully mas- 
sive expression Bf (x, from Eq. ©. In this way, all the mass suppressed effects that 

are not present in Eq. ([1]) but are known from Eq. ([2]) are included. 

In order to do this in a systematic way, we define thus the massless limit of the massive- 
scheme expression Eq. ( [6]) , namely 



L. y \ y m z / 

Jx y i=q,q,g V y 7 

where 

>(°) f, 91 „(»«+i)rn2^ = ( n ^+i) ( o 2 \\ m r!°)'" (, 9. 



(12) 



^ ( *, ^P l+1 \Q 2 )) = E (^ +1) (Q 2 ))" '*f " U ^ ) , (13) 

p=0 



with B^ ,p given by Eq. ()10p . and the sum over p consistently performed including all 

terms up to the order P in a^ l+1 \Q 2 ) to which the massive-scheme expression has been 
determined. 

The FONLL approximation for F 2 is then given by 

F FONLL (x,Q 2 ) = F^(x,Q 2 ) + F^ n '\x,Q 2 ), (14) 
F^(x,Q 2 ) = \f^ +1 \x,Q 2 ) - F^'°\x,Q 2 )] (15) 

where F"' is the massive-scheme expression Eq. ([6]), and the "difference" contribution 
Eq. (|15[) is constructed out of the massless-scheme expression 

i?(n ; +i) Eq. ©, and the 
massless limit F^ 1 ' ^ Eq. ([12]) of the massive-scheme expression. 
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Because of Eq. (jlip . when Q 2 3> m 2 the FONLL expression reduces to the massless- 
scheme one F^+V. When instead Q 2 m 2 the FONLL expressions differ from the 
massive-scheme one through the "difference" term Eq. (|15p . which is however 

subleading in a s (Q 2 ). Note that in the FONLL scheme continuity does not play a role 
in the matching conditions: at NLO continuity ensues accidentally from the matching 
conditions, but at higher orders subleading discontinuities may arise. 

When Q 2 < m 2 , but W 2 > Am 2 i.e. above the threshold for heavy quark production 
there are various options. A simple possibility is to just use the massive scheme result 
Eq. ([6]). It is easy to see that this amounts to replacing the "difference" term in Eq. (1141) 
by 

F^(x,Q 2 ) = [f zmvfn (x,Q 2 ) - @{Q 2 - m 2 )F^°\x,Q 2 )] (16) 

where i ?ZMVFN (x, Q 2 ) denotes the massless calculation with nj = ni + 1 when Q 2 > to 2 , 
and with n/ = n\ when Q 2 < m 2 : this is identical to Eq. (|14p for Q 2 > m 2 , and it reduces 
to the massive calculation for Q 2 < m 2 because the term in square brackets vanishes there. 
Of course, when Q 2 <C m 2 and W 2 <C 4m 2 the massive calculation in turn coincides with 
the massless calculation with rif = ni, so one recovers the same result as the ZMVFN in 
this region. However, if one only wishes to use the result at scales Q 2 which are never 
much lower than m 2 , a simpler option may be to just use Eq. (I14p everywhere. Then, when 
Q 2 < m 2 but W 2 > Am 2 the term in square brackets will contain some small contributions 
in the nt = ni + 1 scheme evolved backwards below the heavy quark threshold. 

The FONLL formula Eq. ()14[) may look similar to a prescription suggested in Ref . [18] , 
and then further discussed in Ref. [24] (in particular Eq. (5) of the latter), sometimes 
referred to as the BMSN prescription, which is also based on the idea of combining com- 
putations performed in schemes which differ in the number of active flavours. However, 
in the BMSN method the issue of using PDF defined in a single factorization scheme in 
all terms is not addressed (unlike in FONLL, where it is accomplished expressing every- 
thing, Eq. ([7]), in terms of a^ 1 ^ and /( n < +1 )). This leads to inconsistent results beyond 
0(a s ), as stated in Ref. [23], where it is argued that the inconsistency is however nu- 
merically small in practice. Also, contributions proportional to the light or heavy quark 
electric charge are not easily separated in the BSMN method, again leading to (possibly 
small [21]) inconsistencies beyond 0(a s ), unlike in FONLL where this separation can be 
treated in a fully consistent way as we will do explicitly to 0(a 2 ) in Sect. 13.21 below. As 
for the comparison of FONLL to the ACOT and TR schemes, we will come back to it in 
the end of Sect. 13.31 after fully specifying the FONLL scheme up to 0(a 2 ). 

Finally, it should be observed that near threshold the F^ term Eq. f)15j> . though 
subleading, could in practice be non-negligible, and it does not provide any information 
because it contains higher-order logarithmic contributions in a region in which L is not 
large. In practice, it may thus be convenient to suppress this term through a suitable 
kinematic factor when Q is near m^, which is allowed without modifying the accuracy of 
the calculation because this term is subleading, as we shall discuss in the next section. 

2.2 Mismatch in accuracy 

In order for Eq. (|14[) and (|15p to work properly, it is necessary to determine F^ ni+1 ^ with 
an accuracy which is at least as high as that used in the computation of F^ ni \ Only in 
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this case, in fact, in the difference expression Eq. (|15p we are subtracting terms that are 
actually present in F^ ni+1 \ 

However, it is possible to generalize Eq. (|15|) in such a way that it can still be used even 
when F^ ni ^ is known with higher accuracy than _F(" i+1 ). For this purpose, it is sufficient 
to retain in F( n ''°) in (|15p only those terms that are also present in F^ TH+1 \ In this case, 
it is no longer true that when Q 2 3> m 2 the FONLL expression reduces to the massless 
scheme one 

F (n l+ 1) 

up to mass suppressed terms only. In fact, F^ 1 ' ^ in (JTSJ) and in 
(|14p no longer cancel in this limit, since some terms have been excluded from F^ Tlh °\ It is 
however still true that the FONLL reduces to the massless scheme one 

F (n l+ 1) 

up to mass 

suppressed terms and terms of higher order in 0(a s ). Further on, we will illustrate an 
application of this approach to the charm structure function, where i^ n ' +1 ) is evaluated 
at NLO, and i^ n! '°) is evaluated at order 0{a 2 ). 

2.3 The heavy quark threshold 

As already stated, when Q ~ m the F^ term becomes totally unreliable, since it contains 
higher-order contributions enhanced by powers of L, in a region in which L is not large. 
It may then be advisable to suppress F^ Eq. f)15|) in the threshold region. Because 
F^ is formally subleading, this of course does not change the nominal accuracy of the 
calculation, but it may in practice improve the perturbative stability and smoothness of 
the result. 

This suppression can be obtained in various ways. Two classes of possibilities which 
have been considered in previous studies consist of introducing a threshold factor or a 
rescaling variable. In the former case, one replaces F^ Eq. f)15|) with 

F^ th \x,Q 2 ) = f thT (x,Q 2 )F^(x,Q 2 ), (17) 

where / t h r is such that F^ d > tfl '(x, Q 2 ) only differs from F^ d \x,Q 2 ) by terms that are 
power-suppressed for large Q 2 , namely 

i 2 ' 

Q 2 



fMx,Q 2 ) = l + 0[^) (18) 



but it enforces vanishing of F^ d,th \x,Q 2 ) below the threshold at Q 2 = m 2 . A suitable 
choice is 

Mr(x,Q 2 ) = ®{Q 2 -m 2 ) fl-V) , (19) 

where the factor in brackets ensure that fthx{x, Q 2 ) and its first derivative with respect to 
Q 2 are continuous at Q 2 = m 2 . 

When using a rescaling variable, one instead replaces in all convolutions which enter 
the expression for the variable x with a new rescaling variable x( x ,Q 2 ), so that the 
factorized expression of all structure functions is now given by 

F^(x,Q 2 )=x f 1 ^c(x^±,a(Q 2 ))f(y,Q 2 ), (20) 

where the rescaling variable x{x, Q 2 ) is chosen in such a way that again F^ d ^\x,Q 2 ) only 
differs from F^ d \x,Q 2 ) by power suppressed terms, namely 
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X ( X ,Q2)= x+ o(^y (21) 

but it is such that, viewed as a function of x, F^ d ' th \x,Q 2 ) only has support above 
threshold. The threshold can be set at the physical production value W 2 = m 2 by choosing 

x = *(i + -Qr). (22) 

This choice was adopted in Ref. [HE], and the ACOT method supplemented by it is called 
the ACOT-x prescription; it is used among others in the recent CTEQ parton fits PQ0 

A generalization of the \ variable was introduced in Ref. [26], by defining a one- 
parameter family of variables which interpolate between x and Xi but all of which set the 
threshold at the physical value. In fact, it was argued in Ref. [15] that a full treatment of 
heavy quark threshold such as ACOT or the FONLL method discussed here could be well- 
approximated by simply using the ZMVFN, but with x replaced by a suitable fine-tuned 
X-like variable. Be that as it may, we note that within the FONLL (or ACOT) method 
there is no conceptual advantage in introducing a suppression at the physical threshold 
W 2 over a suppression at Q 2 = m 2 , given that for most of the interesting kinematic region 
the physical threshold correspond to values of Q 2 which are lower than m 2 , and where 
thus the massive expression is being used anyway. 

In Ref. [17], which dealt with hadroproduction of heavy quarks, yet another form of a 
suppression factor was proposed, based on the transverse momentum of the heavy quark 
pair. We remark here that there is no straightforward generalization of the x-scaling 
prescription to the hadroproduction case. 

All these threshold modifications of higher order terms, being subleading, can neither 
be justified nor be excluded a priori. Their purpose is to mimic in a phenomenological 
way the effect of finite mass effects not included in the calculation. They can thus be 
validated to a certain extent by comparing results obtained using these prescriptions with 
exact results when the latter are known, though of course this does not prove that their 
effect will be the same or similar when the exact result is not known. We will perform 
this validation in section 14.21 below. 



3 Implementation of the FONLL method 

We shall now work out explicitly the FONLL prescription Eq. (fT4l) up to 0(a 2 ). We 
will discuss separately the contributions to the structure function in which the virtual 
photon couples to the heavy quark only ("heavy" structure function Fh) and the rest. 
The basic ingredient in the FONLL construction is the scheme change of Eqs. (|3][3|). We 
will first work this out up to order a 2 , then use it to construct up to the same order the 
various contributions to the FONLL expression of Eq. (|14|) . We will then combine these 

4 Note that a slightly different definition of rescaling is also possible, whereby one simply lets 

p(d,x) =F W(x(x,Q 2 ),Q 2 ). (23) 

This definition, which differs by the above one by a factor J, has been for instance discussed in Ref. [11) ; 
however, the definition Eq. (J20J is used by CTEQ [25]. 
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contributions into expressions for the structure functions. Throughout this section we will 
discuss a generic structure function F(x,Q 2 ), and all formulae will be valid for both F2 
and Fl unless otherwise stated. 

It is important at this point to recall that the meaning of "NNLO" in this context 
is partly a matter of convention: one may use an absolute definition, where LO is 0(1), 
NLO the 0(a s ) and so on, or a relative definition, in which the LO is defined as the first 
non-vanishing order, or a combination of the above according to the observable which is 
being considered. Various options which are relevant for PDF determination up to NNLO 
will be discussed in Sec. 13.31 Finally, we will provide explicit numerical checks of the 
consistency of our procedure in the various cases. 



3.1 Scheme change 

First of all let us discuss the details of the change between the scheme with ni flavours 
(massive scheme) and that with n/ = ni + 1 flavours (massless scheme). For future 
convenience, we define 

a s (Q 2 ) = ^^; /3o = 2tt6o, (24) 

with 

(25) 



33 - 2n f 



12vr 

The running of the coupling with nj flavours to lowest nontrivial order is given by 



O 2 

a s (Q 2 ) = a s (m 2 ) - b log \a 2 s + 0(a\). (26) 



3.1.1 Matching 



At a given value of Q 2 the relation between the coupling constant in the massless and 
massive schemes, Eq. ([3]), reads 

a^ +l \Q 2 ) = a^\Q 2 ) + ^a 2 L + 0{a% (27) 

which immediately implies that the two schemes coincide when the scale is equal to the 
heavy quark mass: 

a^ l+1 \m) = a^\m) + 0{a\). (28) 

Since all coefficient functions at NNLO are at most of order 0(a 2 ), we will not need terms 
of order 0(af ) in the matching condition for a s . 

The relation between parton distributions, Eq. (j4]), for light partons in the two schemes 
at Q 2 = m 2 is also trivial up to order a s , 

f^ +1 \ x ,m 2 ) = f^\ x ,m 2 ) + 0(a 2 ), (29) 

but it becomes nontrivial starting at 0(a 2 ), because at this order contributions from heavy 
quark loops arise, which can thus contribute to the renormalization of light flavours and 
gluons in the massless scheme where n/ = nj + 1. We have thus 

ft l+1 \x,m 2 )=ft'\x,m 2 ) + a 2 J ^ £ K H Wfj^ (f V) + ( 30 ) 
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where the indices i, j take the values q, q, g, and the functions K qq , K gq and K gg were 
computed in Ref. [18], where they are respectively given in Eqs. (B4, B5, B7). On the other 
hand, in the absence of an intrinsic heavy quark contribution, the matching conditions 
Eq. in the case of the heavy quark distribution are 



Jh 



(x,m 2 ) = f^ +1 \x,m 2 ) = 0(aj) 



(31) 



As already mentioned, the coefficient functions C% start at order one only in the case 
of the light quark contribution to F2 (in the absence of intrinsic heavy quark contribu- 
tions). Therefore, up to order a 2 , the only nontrivial contribution to matching Eq. (|30l) 
is due to light quarks: gluon contributions to structure functions start at 0(a s ), thus the 
contribution from K gq and K gg to structure functions only start at O(af). Therefore, in 
this paper we will only use Eq. f)30f) for light quarks, in the form 

if < )(x,m 2 ) = f^ +1 \x,m 2 ) - at f ^K qq (z)f^ g,m 2 ) + 0{a% (32) 

witfH 

K m = C F T R [ilfl (\ log 2 , + § tog , + p) + (i-,)(| tog , + ^ 



(33) 



Eqs. ((30]), is 



3.1.2 Scale dependence 

While the matching condition for parton distributions at the scale Q 2 = m 2 
sufficient in order to fix the relation between PDFs in the massive and massless schemes, 
in order to implement the FONLL approximation we need to express the massive-scheme 
PDFs in terms of the massless-scheme ones at a generic scale Q 2 , as required in Eq. ([6]). 
The relation at any scale can be obtained from the matching condition Eq. (|30p by solving 



the evolution equation for both f^ ni+1 ' and /; in the respective schemes. Because only 
the light quark coefficient function starts at 0{cP s ), up to order a 2 matching conditions at 
0(a s ) will be sufficient for gluons, while 0(a 2 ) will be required for quarks. 
Evolving up to Q 2 Eq. ([29]) to order a s gives 

Jx Z ,: 



(34) 



where Pfj(z) are leading-order Altarelli-Parisi splitting functions in the two schemes, and 
the sum runs over all PDFs. Noting that 



p(ni+l),0, 



(35) 



it immediately follows that the relation between gluon distributions at scale Q 2 is 



& l+1 Hx,Q*) = f^>(x,Q*) 



a s L 2 ^f^(x,Q 2 ) + 0(a 2 ) 



(36) 



5 Note that the expression given here differs by a factor 4 from that obtained setting fi = m in Eq. (B4) 
of Ref. [18] because we expand in powers of ^ and not |^ as in that reference; note also that we have 
absorbed all distributions into the + prescription. 
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To order a 2 Eq. (|34p is modified because the initial condition at Q 2 = m 2 becomes 
Eq. (|30p . and also because the evolution mismatch receives three further contributions: 
from the change in running coupling Eq. (|27p . from the second iteration of the leading- 
order splitting functions P$ , and from the next -to-leading order splitting functions P i j \ 
Combining these contributions, the matching condition Eq. (|30p becomes 

/{"•W) = /^ +1) (x,g 2 )-^|^(z)/^ +1 )g,m 2 ) 

- a ° L f t^i s >< 4 n,+1) (f '^ 2 ) (37) 



^ /lE [^r 1> ' 1 W-- P i"' > ' 1 W] x/j"''(f.<3 2 )+0(aj), 
Z j eggs 



where each line in Eq. (137p corresponds respectively to each of the contributions listed 
above. Note that, because of Eq. (|35|). the 0(a s ) correction vanishes for light quarks. 
Note also that the f g contribution cancels between the second and third line in Eq. (|37p . 
Using the fact that [27] 



Pfr +1) '\z) ~ Pfr^M = -«i<A»(*), (38) 
1 + z 2 /2, 10 \ 4 



with 



A qq (z) = C F T R 



-logz + — + -(l-z) 



(39) 



(m+i) f^ Q2 



l-z\3° 9/ 3 

we reproduce the result of Ref . [18] : 

4 ni \x,Q 2 ) = f^ +1 \x,Q 2 )-a 2 f ^K qq (z,L)f : 

K qq (z,L) = K qq (z) + ^*^p£«>.° W - LA 99 (z), (40) 
with K qq (z) and A qq (z) given respectively by Eqs. 
3.2 The FONLL approximation up to 0(a 



It is now convenient to separate off explicitly the light and heavy contributions to the 
generic structure function F and its corresponding coefficient functions Cf. 

F(x,Q 2 ) = F(x,Q 2 ) + F h (x,Q 2 ); 
Ci(x,a s (Q 2 )) = Cij(x,a s (Q 2 )) + Ci, h (x,a s (Q 2 )) (41) 

where we define Fh and Ci h as the contributions to F and C% respectively which are 
obtained when only the electric charge of the heavy quark is nonzero. Note that therefore, 
up to 0(a 2 ), in the coefficient functions the label I, h denotes the quark to which the virtual 
photon couples, whereas the label i denotes the parton that enters the hard scattering 
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processO The separation of Eq. f)41[) will be carried through also for the corresponding 
FONLL expressions (Eqs. (|14| I15p ). i.e. we will write 

if ONLL (*,Q 2 ) = if%,Q 2 ) + F{ d \x,Q 2 ), (42) 
^ FONLL (x,g 2 ) = Ft\x,Q 2 ) + F^\x,Q 2 ). (43) 

In the following Sects. [3T2. 1113.2.21 we will provide expressions for the various contributions 
to Eq. (HQ up to 0(a 2 ), with more explicit formulae collected in Appendix lAl 

For the rest of this section we shall assume that the generic structure function F(x, Q 2 ) 
refers to the structure function F2(x, Q 2 ) in electromagnetic deep-inelastic scattering. The 
specific differences for the cases of both Fl(x, Q 2 ) and weak-mediated DIS structure func- 
tions will be discussed at the end of Sec. 13.2.21 

3.2.1 Order a s 

At order a s the matching condition for parton distributions Eq. f|30[) is the trivial relation 
fq ni \x,m 2 ) = fq ni+1 \x,m 2 ). This means that at this order all light quark and the 
gluon PDFs should be matched by assuming them to be continuous at Q = m, while 
the (non-intrinsic) heavy quark PDFs should be simply taken to vanish at Q = m. The 
matching condition at a generic scale Q 2 , used to determine the coefficient functions Bi 
of Eq. © in terms of the original massive-scheme ones C- starts at order a s , Eq. (|34[) . 
However, Eqs. (]36( l37"|) show that the correction is nonzero only for the gluon PDF. Since 
the gluon coefficient function starts at 0(a s ), it follows that the coefficient functions Bi 
start differing from the c\ ni ^ only at 0{a 2 )\ 

B t (z, ^ 2 ,4 ni+1 \Q 2 )) = C\ ni) (z, %Mr i+l) (Q 2 )) + 0(a 2 s ). (44) 

At this order all "light" coefficient functions are the same in the massless and massive 
scheme: 

C\f (x,g,a s (Q 2 )) =C^ (x,a s (Q 2 ))+0(a 2 s ) , (45) 

which immediately implies that the FONLL expression of Fi(x,Q 2 ) trivially reduces to 
the massless scheme one: 

F ; FONLL (x,Q 2 ) = Ft d \x,Q 2 ) + F^(x,Q 2 ) = F^ +1) (x,Q 2 ) + 0(a 2 ). (46) 

where we have 

if)(*,Q 2 ) 
F^ l \x,Q 2 ) 

6 Note that up to 0{a 2 a ) the coefficient functions are the sum of contributions which are quadratic in the 
quark electric charge; however at order 0(ot 3 a ), interference terms proportional to the product of different 
electric charges do arise (see Ref. [28]). 



i=h,h 



E / c i7 +1) (-,a s (Q 2 ) ) ft l+1 \y,Q 2 ), 



(47) 



(48) 
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Figure 1: The 0(a s ) diagram for heavy quark production. 



The "heavy" coefficient functions C^h, for i corresponding to any light quark are also 
trivially the same at order 0{a s ) in the massless and massive scheme, since they both 
vanish at that order, but the gluon one is different. We have 

°^ (*' S' as{Q2) ) = *.(<?)*%Cjr ) ' 1 (*, $) + O (cQ , (49) 



where 



r<{ni),l 
^9 



(W 2 - Am 2 ) x T R [(z 2 + (1 - z) 2 + 4ez(l - 3z) - 8e 2 z 2 ) lof 
+(8z(l - z) - 1 - 4ez(l - z))v ] , 



1 + v 
l-v 
(50) 



with 



e = m 2 /Q 2 , v = ^l-4m 2 /W 2 . 



(51) 



The factor of 2 associated with the squared charges in Eq. (149j) accounts for the contribu- 
tion of both the quark and the antiquark. This Feynman diagram for this contribution is 
shown in Fig. [TJ 



The massless limit Eqs. (HUEU) of C$ Eqs. (jUSl is 



B 



(o),i 

9,h 



zX 2 )=2e 2 C^ 



m- 



where, according to Eq. (|10p the index 1 denotes the first order in a s , and 



{z* + (l-zY)]o, 



2N1 _Q 2 (l-z) 



m 2 z 



+ (8z(l-z) - 1) 



(52) 



(53) 



which in the limit Q 2 = m 2 reproduces the massless scheme coefficient function: 

clX 1 \z,a s (Q 2 ))=a s 2elc^(z,l) . 



(54) 
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The FONLL expression for the "heavy" component is given by 

if ° NLL M 2 ) = if* W) + ifW)- (55) 



The two contributions on the right-hand side of Eq. (|55p are respectively the massive- 
scheme heavy-quark structure function, given by 

Ft\ X ,Q 2 ) = xf^C^ ( X -S a ^) /i^W), (56) 

with Qp,a s ^J given by Eqs. (|49l50j) : and the "difference" contribution 

Fl d \x,Q 2 ) = x f d ^\c { X l) (-^s{Q 2 )) [ft l+1) (y,Q 2 ) + t l+1 \y,Q' 



(n,+l) / X 



C g,h 



y 



« s (Q 2 ))-<l(^' as(g2) ))^ +1)(y ' Q2) 



(57) 



where C^" 1 " 1 ^ and B^ h are given by Eqs. (f5~ij) . ([52]) . respectively, and C^ +1 ^ is the 0(a s ) 
quark coefficient function for production of a (massless) quark of charge e^. 

It is easy to see explicitly that, in the region where L is not large, the "difference" term 
is of order 0(a 2 ), so that to 0(a s ) the FONLL expression coincides with the massive- 
scheme one. Indeed, Eq. (|3ip implies that the heavy quark distribution is 0(a s ). Use of 
the leading-order QCD evolution equations immediately leads to 



h(y, Q 2 ) = h(y, Q 2 ) = ^f±L J ^T f (z 2 + (1 - z) 2 )g Q 2 ) + 0(a 2 s 



(58) 



Noting that C Q: h {x, a s (Q 2 )) = 5(1 — x) + 0(a s ) it immediately follows that the two terms 
in Eq. ([57|) cancel up to terms of order 0(a 2 ). 

3.2.2 Order a 2 s 

We will now work out the FONLL expressions for structure functions at 0(a 2 ). We expand 
structure functions and coefficient functions in powers of ai n ' +1) (Q 2 ), Eq. PU, as follows: 

p /J ni+1) ((l 2 )V P ( ry {nt+1) (<l 2 )Y 

p=o V 27r / p=o V 27r / 

and similarly for other quantities. In order to deal with manageable expressions we will 
collect in Appendix [A] the expressions for the 0(a 2 ) coefficients Cf, f f °nll,2^ an( ^ go Qn ^ 
while the contributions up to 0(a s ) were given explicitly in Sec. 13.2.11 



At the NNLO order, both "light" and "heavy" coefficient functions differ in the massive 
and massless scheme. As far as "light" coefficient functions are concerned, the gluon one 
is the same 

m 2 
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but the quark one O-^g'i is different, because of the contribution coming from the emission 
of a gluon which in turns radiates a h—h pair, shown in Fig. [2j 

Further contributions may arise when re-expressing the strong coupling and parton 
distribution in terms of their massless-scheme counterparts for quarks we find that the 



0{a%) contribution Bf l for i / g (Eq. ©) is given in terms of the coefficients (Eq. (|59 
of the expansion of the massive-scheme coefficient function by 



where C^' 2 is the coefficient function for the i— th light quark or antiquark, K qq (z,L) 



(ni),2 



m 1 



e*K qq (z,L) 



2T R 



LC, 



(ni),l 



(61) 



is given in Eq. (|4(jp and cf 11 '' 1 is the standard [27j (massless) NLO contribution to the 
coefficient function for the i— th quark or antiquark. The term proportional to C\ 

(using Eq. ([27]) ) in the 0(a s ) correction 
for i ^ g is given by 



(m) 



m terms of a s 



arises when expressing a 
to C\ ni \z). The massless limit of Bf h Eq. (JBT 



B 



m- 



C 



(n f ,0),2 
i, / 



eiK w {z,L) + -^LC t 



1T R ( ni ),i 



(62) 



where C; j 2 denotes the massless limit of the coefficient function in the massive scheme. 

In the gluon sector, even though the coefficient function C„l coincides with its 
massless counterpart, according to Eq. (|60p . both the expression in terms of their massless- 
scheme counterparts of the strong coupling Eq. (j27[) and of the gluon distribution Eq. (1361) 
contribute to the determination of Bg (Eq. Q), as discussed in Sec. 13.1.21 The former 
contribution is of course the same as in Eq. (|6ip . and the latter turns out to cancel it 
exactly, so that in the end we get 



B 



9,1 



Q 1 



777/ 



(63) 
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Figure 3: The light-quark initiated 0(a 2 s ) diagram for heavy quark production which contributes 
to the "heavy" coefficient function Fh Eq. (|4Tj). 

namely, the coefficient function is the same as in the massless scheme after all. 

Turning now to "heavy" coefficient functions, both the quark and gluon ones, C^^' 2 , 
with i = q, g are nontrivial in the massive scheme, because heavy quarks contribute both 
to quark-initiated (Fig. [3]) and gluon-initiated processes (Fig. H]). When re-expressing in 
terms of the massless-scheme strong coupling and PDFs, we now find 



Indeed, for an incoming light quark the "heavy" contribution only starts at 0(a 2 ), so up 
to 0(a 2 ) it is unaffected by changes in the coupling and PDFs, while for an incoming 
gluon the argument leading to Eq. ([63]) also applies, and leads to the same conclusion. 

Collecting this information, we can construct the FONLL expression for the various 
contributions to the structure function. For the "light" component Fi(x, Q 2 ), we find that 
the 0(a 2 ) contribution to the FONLL expression is given by 




(64) 



F, 



,FONLL, 2 



(x,Q 2 ) = F^ 2 (x,Q 2 ) + F^\x,Q 2 ) 



(65) 



Here the contribution to the massive-scheme structure function is 




(66) 



where the sum runs over light quark flavours and antiflavours and Bf, and /j are 
respectively the massive coefficient function Eq. (|6ip and PDF for the i— th light flavour. 
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The 0(a 2 ) term of the "difference" contribution is then given by 



*f (*.<J S 



z E 

ij^h,h,g 
i=h,h 



1 dy 

y 



c 



(rn+l),2 ( X 



i: I 



D 



(0),2 fx 



ij 



y m' 



f, 



(ni+l) 



(y,Q 2 ) 



dy_ c (ni+l),2 f x\ f ( ni+1 ) 



i. I 



fr^>(y,Q 2 



(67) 



where B^' 2 are given by Eqs. <\62\ [63]) . and the NNLO massless light quark and gluon 
coefficient functions C , ^ +1 ' 1 ' 2 , first computed in Refs. [29j[30] are respectively given in 



K+l),2 
i, h ' 

Eqs. (B.2-B.4) and Eq. (B.6) of Ref . [31] HI 

Now that the structure function is computed to 0(a 2 ), when L is not large, the 
difference contribution F d must start at the order a^. The second term in Eq. (|67p 
manifestly starts at this order, since it is proportional to the product of the heavy flavour 
parton density, which starts at order a s , and the coefficient function for an incoming 
heavy quark to produce a light quark that couples to the photon, which is of order a 2 . In 
order for the first term in Eq. (j67[) not to be of 0(a 2 ), the difference of coefficients in the 
square bracket must vanish identically. Notice therefore that in actual fact the coefficient 
^f 1 } 2 (§ ' ^) does not depend on Q 2 /m 2 . In order to verify this cancellation, we rewrite 
the first term on the r.h.s of Eq. (|62p as 



C 



(ni,0),2 



y(ni),2 l 



) + e 2 C^ 2 (zX 2 



(68) 



where C^ l ^' 2 (z) stands for the MS coefficient function of flavour i with n/ light flavours, 



and 



r NS,2 f Q 2 \ _ 1 [ r NS,(2) f ' m2 \ , r 



NS,S+V,(2) 

q 



m 

Q 2 " 



(69) 



7 Note that in Ref. [31] coefficient functions are given as a series in powers of j£ 
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Figure 5: The light-quark initiated Q(oq) diagram in which heavy quarks contribute to the gluon 
self-energy. 



where L^ s is given in equations (D.8) and (D.10) of Ref. [32] and the factor of 1/4 accounts 
for the fact that ol s /(Ait) rather than a s /(2ir) is used in that reference. In other words, 
Eq. (|68p splits the coefficient functions into terms not involving heavy flavours at all, and 
terms involving a heavy flavourjl like the diagrams shown in Figs. I2|5I 

Since the C^ l+1 ^' 2 coefficients are also MS coefficients (and since they depend linearly 
upon n;), we have 

C% l+1) '\z) ~ C\f>\z) = ±-C%' +1) > 2 (z), (70) 

so that, in order for the first term of Eq. (|67|) to be of order we must have 

_d_ 

dni 



-Ci n > +1) > 2 (z)-e 2 C NS > 2 (z,m 2 /Q 2 ) + 



'4K qq { Zl L) + 2 ^LC\ n ^\z) 



0. (71) 



Using Eqs. gQD, (D.8) and (D.10) of Ref. [32], and Eq. (4) of Ref. (29] for the coefficient 
of ni in C^™ !+1 ^' 2 , we can easily verify Eq. ([711) 1^1 

The structure of the C(a 2 ) contribution to the FONLL expression for the "heavy" 
component of F is similar to that of its 0(a s ) counterpart Eq. (155p . namely 



F™' 2 (x,Q 2 ) = Ff' 2 (x,Q 2 ) + F^^ 2 (x,Q 2 ), (72) 
except that now both light quarks and gluons contribute to the massive-scheme term: 



Ft l) '\x,Q 2 ) = x f 1 ^ £ B ^' 2 ft l+1) M 2 ), (73) 

Jx y .^fr ' \y m 2 J 



8 By "terms not involving heavy flavours" we also mean here terms where the heavy flavour contribution 
totally decouples in the decoupling scheme, which is typically the case of heavy flavour loop corrections to 
on-shell gluon propagators. 

9 We remark here that the first moments of C7 ( "< +1) ' 2 , C§ l) ' 2 , C NS ' 2 (z,m 2 /Q 2 ), K qq (z,L) and c\ ni)A 
all vanish individually. This fact, which can be explicitly verified, is a consequence of the Adler sum rule. 
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with the sum running over all light quark and antiquarks and gluons, and B^J 1 ^' 2 given by 
Eq. (I64p , The real diagrams corresponding to the gluon initiated contributions to Eq. (1731) 
are shown in Fig. HI while the corresponding light-quark initiated contribution is shown 
in Fig. El 

The difference contribution is now given by 

{d\2, n2 x _ „ ST f 1 d yMni+l),2 ( x\ (n l+ \) 



E 

i^h,h,g 



'dy 
y 



\y J \y m 



ft l+1 \y,Q 2 )- (74) 



Here C^ +1) ' 2 are the massless-scheme 0(a 2 ) contributions to the standard coefficient 
functions for production of a quark of electric charge eh, from Refs. [29H31] and already 
used in Eq. (|67p . with the heavy quark treated as another massless flavour, while -B^' 2 
are the massless limits of the massive coefficient functions Eq. (|64p . given in Eqs. (D4,D6) 
of Ref. [32] for gluons and quarks respectively. 

As mentioned in the beginning of Sec. 13.21 all results so far apply to the structure func- 
tion Fi(x, Q 2 ). However, they also hold for the longitudinal structure function Fl(x, Q 2 ), 
whose coefficients are given in the references cited above along with those for Fi- The 
only difference is that, because the structure function Fl(x,Q 2 ) starts at 0(a s ) also in 
the light quark sector, for this structure function one finds the simplified expression 

BhM $) = C»* U $) - ^L C f>-%) (75) 



instead of the more complicated relation Eq. (1611) . 

While we have provided up to now explicit results for electromagnetic DIS structure 
functions, it is straightforward to generalize the above discussion to the case of weak- 
mediated neutral current DIS: indeed, the only modifications consists of the replacement 
of the electromagnetic charges by the corresponding weak charges. 

The generalization to charged-current (CC) DIS is slightly more complicated. Firstly, 
full 0(a 2 ) CC coefficient coefficient functions are not available: only their asymptotic 
Q 2 3> m 2 limits [33] and their threshold behaviour [34] are currently known, so that the 
FONLL can only be implemented to this order. For CC scattering initiated by a heavy 
quark, like strangeness production though W— boson exchange, the FONLL formulae of 
Sec. 13.2.11 straightforwardly generalize, taking into account only the different charges and 
kinematics. On the other hand, for CC scattering with a light quark in the initial state, 
such as charm production in neutrino scattering (commonly known as "dimuon produc- 
tion") the heavy quark PDF does not enter, so the FONLL results reduces to the massive 
calculation (with scheme ambiguities only entering at 0(a 2 )). 



3.3 Structure functions: schemes and perturbative ordering 

Various possibilities for the definition of perturbative ordering are possible and have been 
advocated, especially in the context of parton fits including DIS data. This is due to the 
fact that LO parton distributions can be used for the computation of any hard process 
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at the first nonvanishing oder: it is then natural to adopt a "relative" definition of per- 
turbative ordering, where LO is the lowest nontrivial order at which the process starts 
occurring. This is to be contrasted with an "absolute" definition of perturbative ordering, 
such that LO refers to 0(a®) (parton model), NLO to 0(a s ), and so forth. The issue 
is nontrivial because in a global parton fit one may want to combine data for the total 
structure function F, which starts at 0(a®), and those for the heavy structure function 
Fh, which starts at 0{a s ). 

One option is then to simply adopt an absolute perturbative ordering, whereby at NLO 
all quantities are computed at 0(pq). This choice has been adopted for instance in the 
CTEQ parton fits [T]. We will refer to it as scheme A henceforth. However, as discussed 
in Sec. I2.2l we can also consistently use different accuracies in the computation of different 
contributions. At NLO we can then compute the total structure function at 0(a), but 
use the 0(a 2 ) results of Sec. I3.2.2l for the determination of Fh (both with NLO PDFs). In 
such case, it is not necessary to compute any of the massless-scheme coefficient functions 
C\ ni+l) beyond NLO (i.e. beyond 0{a s )) because this does not improve the accuracy of 
the NLO calculation in the massless sector, but all massive coefficient functions B 2 and 
their massless limits B^' 2 should be computed to 0(a 2 ) in order to have this accuracy 
in the massive sector. In order to ensure consistency of the subtraction it is sufficient 
to include only logarithmic terms in B^' 2 , because the corresponding non- logarithmic 
terms are not included in the massless coefficient functions cf ni+ . We will refer to this 
as scheme B. 

This scheme B is reminiscent of the method proposed in Refs . [HHlO|. and used e.g. 
in MSTW PDF fits [14] : there too, the 0(a 2 ) expression of the massive cross section is 
used in conjunction with NLO parton densities and coefficient functions, although within 
a totally different approach. 

We can then also pursue a full 0{a 2 s ) computation, by simply performing the absolute- 
orderd computation of scheme A to one extra order, with NNLO PDFs. We will refer to 
this as scheme C. However, we cannot pursue scheme B to one extra order because O(a^) 
massive coefficient functions are not known yet, though their large Q 2 limit was recently 
determined in Ref. |35j . In NNLO MSTW fits |14j they have been modelled based on their 
known large- and small x behaviour. 

In summary, we will consider three options for perturbative ordering, whose consistency 
will be checked explicitly in the following Sec. 13. 4| and whose phenomenological implica- 
tions will be discussed in Sec. HI namely: 

• Scheme A: the 0{a s ) FONLL expressions of Sec. 13.2.11 with NLO parton densities. 

• Scheme B: the 0(a 2 ) FONLL expressions of Sec. 13.2.2] but using NLO parton densi- 
ties and 0(a s ) massless coefficient functions, and retaining the 0(a 2 ) contributions 
to the massive coefficient functions. In this case, the massive expression exceeds in 
accuracy the massless one when L is not large. Thus, according to section 12.21 in 
the massless limit expression of Eq. (|15p . the non-logaritmic 0(a 2 ) term in F^ n '' ^ 
should not be included. 
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Q/m Q/m 



Figure 6: The F2 c structure function computed in the massless scheme Eq. [T] (ZM) and in the 
massless limit of the massive scheme Eq. (|T^|) (MO), in the three perturbative ordering schemes 
A-C of Sec. 13.31 The structure function is plotted as a function of Q/m, for x = lO - ^" -1 ^, with 
integer 1 < j < 5 (from top to bottom at large Q/m). The plots on the right show a magnification 
of the threshold region. 
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Figure 8: The Fi c structure function computed in the massive scheme Eq.[2](M), masslcss scheme 
Eq.[T](ZM), massless limit of the massive scheme Eq. (JT2J (MO) and FONLL scheme Eq. ([H)>. The 
schemes adopted in the three columns correspond to the different combinations of perturbativc 
orders discussed in Sec. 13.31 Results at the scales Q 2 = m 2 . = 2 GcV 2 and Q 2 — 4 GeV 2 are given 
in the two rows. 



• Scheme C: the 0(a 2 ) FONLL expressions of Sec. 13.2.21 with NNLO parton densities. 

We address now the question of the relation of the FONLL scheme with other ap- 
proaches. As already said, our B scheme is reminiscent of the TR scheme, although it is 
certainly not identical to it. The complexity of the TR implementation has prevented us 
from obtaining a clearer view of the differences. 



Our A scheme is easily shown to be equivalent to the SACOT variant of the ACOT 
scheme. If a x- scann g prescription is also applied (with the same definition used by the 
ACOT group), the A scheme becomes equivalent to the SACOT-x one. There is no analog 
of our B scheme within ACOT. A full discussion of ACOT at NNLO as not yet appeared 
in the literature. We believe, however, that our C scheme should be equivalent to SACOT 
at the NNLO level, or at least we do not see any reason why it should differ from it. A 
benchmarking of the FONLL, ACOT and TR prescription, which in particular confirms 
that FONLL and SACOT-x agree to NLO, has been recently performed in Ref. |36j 
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Figure 9: Same as Fig. [51 but at the higher scales Q 2 = 10 and Q 2 = 100 GeV 2 . Note that a log 
scale is now adopted in the y-axis. 

Recently, the BMSN prescription which we briefly discussed in Sect. 12. ll after Eq. (|16D . 
has also been used for PDF determination in Ref. [19]. On top of the differences already 
mentioned there between the FONLL and BMSN approaches, in Ref. [19] the resummation 
of powers of L Eq. ([5]) which we discussed at length in Sect. 12.11 is not performed [20J: 
this may be adequate in PDF determinations which do not use data at large Q 2 , but it 
requires a subsequent transformation to a scheme where such a resummation is included 
for the sake of phenomenology at large Q 2 , as indeed done in Ref. [19] . 



3.4 Consistency of the FONLL procedure 

The consistency of the FONLL method is guaranteed by the subtraction of the overlapping 
terms from the massless result, as implemented in Eq. (|14p . Ideally, one should be able 
to verify that, in the m — > limit and with light flavour PDFs fixed at the scale Q 2 , the 
massless scheme (ZM, henceforth) result Eq. ([1]) and the zero-mass limit of the massive 
(MO, henceforth) result Eq. (|12[) should coincide up to higher orders in a s . This is in 
practice quite difficult in general [37] ; however, several features that must hold in the 
relation between the ZM and MO result near threshold are rather easier to verify: in 
this section we will perform some of these checks, thereby verifying the consistency of 
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Figure 10: Same as Fig. HI but for the structure function Fl c - 



the FONLL procedure. The discussion which follows applies to a generic deep-inelastic 
"heavy" structure function Fh(x,Q 2 ) Eq. (|4ip : we will show explicit results for both the 
structure functions F 2 h(x,Q 2 ) and Fi^x^Q 2 ). 

Recall how the perturbative expansion and the logarithms of Q 2 /m 2 are organized in 
the ZM and MO result. The "heavy" structure function Fh(x,Q 2 ) in the massless scheme 
Eqs. (Ill4ip is given by 

Ft' +1 \x,Q 2 ) = x f d -^\c^ l) (-,a s (Q 2 )) (h(y,Q 2 ) + h(y,Q 2 )) (76) 

+ C{ X 1] ^» s (Q 2 )^( q (y,Q 2 ) + q(y,Q 2 )) +C^ +1) (~^(Q 2 )) 9(y,Q 2 ) 
where the choice of scheme for a, and PDFs is immaterial. 



For the sake of a comparison between the ZM and MO results, the heavy quark PDFs 
h,h can be viewed as a function of L = log (Q 2 /m 2 ) , with the light parton densities 
9(2/) Q 2 ) an d q(y,Q 2 ) treated as constants. At NLO we get 

h(y,Q 2 )^ hf\y)a\U+ £ ^(j,)^" 1 . (77) 

i=l, 00 i=2,oo 
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Decomposing the coefficient functions as 

n {ni+l) _ n {ni+l)fl n ( ni +l),l n {ni+l) _ n [ni+l),\ n (ni+l) _ n , 7H \ 

L h,h ~ C h,h + a * C h,h > C g,h ~ a s L g,h ' C q,h ~ U ■ V ' °J 

the ZM structure function is, up to NLO 

j x y y \y / 

+ c tl +1),X (~) {h ( ?(y)<x 2 s (Q 2 ) + 4%)^(Q 2 )) l + c^ +l) ' 1 (~) a.W'Mv, Q 2 ) + o (l 2 ) 

On the other hand, the ZM result can be written as a power series in L by expanding 
the coefficient function as 



(79) 



N 



C^'°) (x,L,a,(Q 2 )) = J2Y, C{ni '° ),i,j W<Q 2 Y Lj - 

i=l j=0 



(80) 



In scheme A of Sec. I3.3| which corresponds to using NLO with fixed-order results 
consistently included up to 0(a s ) result, one gets 



+U g,h 



(n,,0),l,l / X 



a s (Q 2 )g(y,Q 2 )L 



(81) 



Therefore, because of Eq. ([54]) the ZM result Eq. flZ9J) and the MO result Eq. (jHI]) must 
coincide when L = 0, while their slope in L near L = should differ by terms of order a 2 , 
which are present in the ZM Eq. (I79p but not in the MO result Eq. (]8ip . This is borne 
out by the plots Fig. [6] for F2 C (x,Q 2 ) and Fig. [7] for Fl c (x,Q 2 ), in which the structure 
function is plotted as a function of Q/m for several values of x (with m = m c = y/2 GeV). 
We have also checked that the difference in slope is indeed of order a 2 , by repeating the 
plot with different values of a 2 . 

In scheme B, which corresponds to NLO but with 0(a 2 ) contributions included in 
massive coefficient functions, the MO result becomes 



F 



(n u 0) 



(x, Q 2 ) = x 



dy 
V 



c iS' 0),1, ° (f ) ^Q 2 )g(y,Q 2 ) 



+C. 



(ni,0),l,l / X 



g,h 



-) as (Q z )Lg(y,Q z ) + C i 



(n,,0),2,l / X 



g,h 



-)a s (Q 2 ) 2 Lg(y,Q 2 " 

v 



+c 



(m,o),2,i / x 



± a s (Q 2 ) 2 L £ (q(y, Q 2 ) + q(y, Q 2 )) + O (L 2 ) 



(82) 



which now also receives a contribution from the light-quark initiated terms of Fig. [21 which 
start at O (a 2 ) . Note that, as already mentioned in Sec. 13.31 in scheme B non-logarithmic 
O (ot 2 s) are not included in the MO result, because they must match the ZM result, which in 
this scheme is determined at NLO, so that beyond O (a s ) it only contains logarithmically 
enhanced terms. Comparing the MO result Eq. (182p to the ZM one of Eq. (|79|) one sees 
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Figure 11: Same as Fig- El but for the structure function Fl c . 



that in scheme B not only the value at threshold but also the slope should agree, up to a 
small correction O (a^) due to contribution proportional to in Eq. (|79[) . This is borne 
out by the plots in Figs. 

The same reasoning can be pursued to NNLO. In this case, we only have results in 
scheme C, in which fixed-order results are included up to 0{a 2 s ). In this case the ZM and 
MO results differ at threshold both in value and slope through terms of order arising 
both from the PDF and the coefficient function. This mismatch is apparent in Figs. ErEl 
and turns out to be quantitatively different in F2 C and in Fl jC . Indeed, for F2 C (Fig. [6]) 
the difference in slope is very visible, but the difference in value at L = is very small 
and difficult to appreciate. On the other hand for F^ c (Fig. [7]), higher order corrections 
are larger and the differences also in values at threshold are more visible. 



4 Phenomenological impact of the FONLL method 

We will now examine the phenomenological impact of the FONLL method on structure 
functions. We will consider the case, both for i*2 and Fl, of the total structure function, 
as well as the "heavy" contribution Fh to structure functions in the particular case of 
charm, namely F2 C and Fl c - The purpose of our study is to assess the impact of heavy 
quark terms at various perturbative orders, to compare and validate different prescriptions 
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for the treatment of sub leading terms close to threshold, and finally to provide a general 
assessment of the impact of heavy quarks on structure functions. A study of the impact 
of heavy quark terms on a determination of PDFs should be part of a global PDF fit and 
it is thus beyond the scope of our work. 

In order to ease comparison with other studies, we will use the PDFs of the standard 
reference Les Houches [38] set at the scale Qq = 2 GeV 2 , evolved using the QCDNUM [39J 
package, and tabulated in Refs. [381140] both in the massive and massless scheme (with, 
in the latter case, PDFs matched according to the results given in Sec. I3.1.ip . Massive 
0(a 2 ) contributions to coefficient functions have been determined using the numerical 
implementation of the results of Ref. [21], and their massless limit using the numerical 
implementation of the results of Ref. [32] . All the analytic expressions which we have used 
and the computer codes used to implement them are illustrated in the Appendix [A] As in 
Ref. [38], we take the charm mass at m 2 = 2 GeV 2 and a s (m 2 ) = 0.35, which corresponds 
to values in the ballpark a s (M 2 ) ~ 0.120 at NLO and a s (M z ) ~ 0.110 at NNLO according 
to the scheme used (see Refs. [381HQ] for the exact values). We will furthermore neglect 
the b threshold, by letting mi, — > oo. 

4.1 Comparison between FONLL, massive and massless results 

At first, we compare structure functions in the massive scheme Eq. ([2]) (M, henceforth), in 
the massless scheme Eq. ([I]) (ZM, henceforth), in the massless limit of the massive scheme 
Eq. (|12p (M0, henceforth), and using the FONLL prescription Eq. (|14p . We begin by 
showing in the first three plots of Fig. [8] the structure function F2 C , computed using the 
three different combinations of perturbative orders discussed in Sec. 13.31 and referred to 
as schemes A-C, at the charm mass scale Q 2 = m 2 = 2 GeV 2 . We see that in all cases 
the ZM and M0 results either coincide, or else are very close, as we have already seen in 
Sec. 13. 4[ It follows that the "difference" contribution Eq. (|15j) is vanishingly small, and 
thus FONLL result coincides with the massive calculation. 

However, the same quantities at Q 2 = 4 GeV 2 (also displayed in Fig. [8|) show remark- 
able differences. In the scheme A, the ZM and M0 results differ considerably, so FONLL no 
longer simply reduces to the massive result, and it is in fact rather closer to the ZM result. 
This is not the case in scheme B, in which the ZM and M0 results again nearly coincide. 
The case of scheme C is similar to that of scheme B for x > 10~ 3 , but for smaller values 
of x the ZM and M0 result again differ considerably. In Fig. [9] we see that this pattern 
persists at higher energies, although with a reduced impact of mass effects. This means 
that inclusion of the 0(a 2 ) contribution to the massive terms in scheme B significantly 
reduces the impact of the unreliable 0(a 2 ) terms in the difference contribution which are 
present in scheme A, and whose effect persists even at larger scales. However, in scheme 
C further 0{a^ s ) terms are introduced, which may become large at small x due to large 
small x logs [4T] . 

From the plot on the first column of Fig. [U we also see that the difference between 
the M0 and M curves (i.e. the impact of mass effects) in the first order heavy flavour 
coefficient is quite large even at a scale of 10 GeV 2 . This difference seems to be reduced 
in the third column, which suggests a cancellation of mass effects between the first and 
second order in a s . This cancellation should be considered accidental, and it cannot be 
taken as evidence that mass effects are negligible at this scale. 
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Figure 12: Same as Figs. MM but for the total structure function F 2 . 
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We now turn to the structure function Fl c - Also in this case, right at threshold 
(Fig. [T0|) the ZM and MO results compensate exactly in the A and B schemes, though now 
in scheme C they start to deviate at x < 10~ 2 , and only coincide for larger x values. As the 
scale is raised to Q 2 = 4 GeV 2 , the ZM and MO results start to differ, though by a much 
smaller amount than in the case of Fi c . In this case, however, the ZM and MO results at 
and close to threshold are much larger in magnitude than the M result. This means that 
for Fl jC the massless approximation is completely inadequate close to threshold. This is 
clarly seen in Fig. [10\ where the massless result at threshold in all cases A-C is totally 
different from the massive one, to which instead FONLL essentially reduces in cases A-B, 
because of the complete cancellation between the M and MO terms. In case C, however, 
the M and MO terms do not cancel completely at threshold (recall Fig. [JJ especially the 
lower right plot) and thus, because of the the large relative size of these (subleading and 
unphysical) terms in comparison to the massive result, the FONLL result ends up having 
also an unphysical behaviour at threshold. This means that in this case it is crucial to 
implement a suppression of these subleading terms as discussed in Sect. 12.31 and that this 
suppression must be very strong at threshold (as it is using the threshold factor Eq. (I19j) ): 
only in such case the difference contribution Eq. f)15|) is removed and FONLL reproduces 
the M result again. Obviously, this means that the impact of the threshold prescription 
for Fl in scheme C is very strong. Once again, at Q 2 = 10 GeV 2 (Fig. [TT]) the M and MO 
results differ considerably. Even at Q 2 = 100 GeV 2 mass effects are not at all negligible 



In Figs. [T2113I we turn to the to the total structure functions, i.e. the sum of the light 
and heavy contributions Eq. (|4ip . As discussed in Sec. I3.2[ the FONLL prescription has 
no effect on F\ to 0(a s ), and thus is scheme A. Its effect to 0(a 2 ), which comes from 
the Feynman diagrams of Figs. [2][5] is very small: the contribution from these diagrams 
is small in the first place, and its contribution to the FONLL prescription is further 
suppressed by the cancellation Eqs. (|70H71[) . Therefore, the dominant effect of the FONLL 
prescription is on Fh, and thus diluted in the total structure functions. Nevertheless, the 
main qualitative features remain the same. Namely, in scheme A a significant difference 
between the massive and FONLL results arises just above the threshold, and persists up 
to scales as large as 10-20 GeV 2 , with the FONLL result rather closer to the massless one 
for F<l and to the massive one for Fl. This difference is driven by subleading (unreliable) 
threshold contributions, and thus sensitive to the threshold treatment discussed in Sec. 12.31 
to be studied phenomenologically in Sec. 14.21 Its effects are greatly reduced in schemes 
B-C (but with some small x instability in scheme C) and in general for the structure 
function Fl- 

4.2 Comparison and impact of threshold prescriptions 

As explained in Sec. 12.31 the purpose of threshold prescriptions is to modify subleading 
contributions to the FONLL expression which are certainly not accurate but may be 
non-negligible close to threshold. As such, from a theoretical point of view all threshold 
prescriptions (including the choice of not implementing any threshold prescription) are 



for F Lc . 
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Figure 13: Same as Figs. [TUHTT1 but for the total structure function F^. 
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0.20 




Figure 14: The 0(a s ) contribution to the charm structure function, F 2c (x, Q 2 ) with m 2 = 2 GeV 2 
plotted as a function of Q 2 at fixed x at small x. The four curves correspond to the the full massive- 
scheme (M) result Eq. ([BJ , its m — > limit (MO) Eq. ( [12]) , and the results obtained applying to 
the MO a damping factor Eq. (TT5)) or x~scahng Eq. (|2U|) . 

equally justified, and varying the threshold prescription merely provides an estimate of 
unknown subleading terms. However, when the higher order exact result is known we may 
validate the threshold prescription by comparing the subleading terms which are affected 
by the threshold prescription with their exact form. Here we perform such a validation; 
then we assess the impact of different choices for the threshold prescription on the final 
FONLL result. 

The terms which are affected by the threshold prescription are the difference between 
the ZM and MO expressions: they are nonzero because the ZM contains contributions to all 
orders, while the MO is determined at fixed perturbative order. In the simplest (trivial) case 
in which the structure function is determined at 0(a^), the massive contribution to the 
FONLL expression vanishes: the FONLL expression then reduces to the massless scheme 
one which we know to be inaccurate close to threshold. A threshold prescription applied 
at this order might then reduce the difference between the massless scheme expression and 
one in which mass effects are properly treated. In fact, in Ref. [15] the use of a massless 
scheme supplemented by the x _ scaling threshold prescription was advocated as a possible 
approximation to schemes were mass effects are fully included such as the FONLL scheme 
discussed here. 

Of course, the first massive contribution appears at 0(a s ), so in order to assess whether 
indeed a threshold prescription does help we should compare the exact 0(a s ) massive (M) 
result, i.e. F^ ni \x,Q 2 ) Eq. © to its massless (MO) approximation F ( - nu °\x, Q 2 ), which 
is included in the massless (ZM) result, and to the putatively improved versions of the MO 
which are obtained by supplementing it with a threshold prescription, such as the damping 
factor Eq. ([18]) or x~scaling Eq. ()20|) . This comparison is performed in Figs. [T4ll5l for the 
i<2c structure function. We see that at very small x mass effects are small, but that for 
x ^ 0-01 close to threshold the deviation between the M result and its MO approximation 
becomes significant, the full result being suppressed by mass effects in comparison to its 
massless approximation. The threshold factor reproduces well this suppression especially 
at smaller x values, while x~ sc aling provides a bit too much suppression and accordingly 
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Figure 15: Same as Fig. [TJ] but at large x. 



a somewhat worse approximation. That x~scaling provides an excess of suppression for 
i<2 was already noticed in Ref. [26], where modified, improved x variable were suggested 
to remedy this situation. It is interesting to observe that the slightly different x _sca hng 
prescription Eq. (f23|) (discussed in Ref. [2]) leads instead to rather less suppression, and 
in particular, at small x, less than the damping factor. 

Our results so far, Figs. [T41151 show that indeed a threshold factor may help, but 
they are of academic interest as soon as we implement the FONLL prescription at lowest 
nontrivial order, because 0(ot s ) terms are then treated exactly. We thus turn to the more 
practically relevant case of the 0(a 2 s ) terms: these are treated approximately if we adopt 
the NLO FONLL method of scheme A. In Figs. [Toll 171 we repeat the comparison of the M, 
MO, and threshold-corrected MO results, but now for the 0(a 2 ) contributions to i*2 C - Now 
the difference between the exact M and approximate MO result is significant even at very 
small x. It is still true that, while results with threshold prescriptions are closer to the 
exact one, x _sca hng provides a somewhat excessive suppression. However, now at large x 
the damping factor provides insufficient suppression, so the quality of the approximation of 
either threshold prescription is similar, though perhaps the damping factor is still slightly 
better on average. 

If the FONLL method is treated including massive contributions up to 0(a 2 ), i.e. 
using schemes B or C the threshold prescription only starts affecting O(o%) terms. We 
can then no longer assess its accuracy, because exact O(a^) massive coefficient functions 
are unknown. 

We can now study the impact of various threshold prescriptions on the FONLL result. 
This is done in Figs. [T81I191 where we compare respectively the FONLL expression for the 
structure functions F^ c and Fl c , determined at different perturbative orders, i.e. within 
schemes A-C, and either without any threshold prescription, or with either of the two 
threshold prescriptions that we discussed so far. The effect of the threshold treatment is 
most visible just above threshold at Q 2 = 4 GeV 2 , as expected based on the discussion of 
Sec. 14.11 and Figs. [9][8j At this scale, in scheme A the effect of the threshold suppression 
procedure is visible for all values of x; in schemes B-C, where the threshold treatment only 
affects O(a^) terms, the effect is generally small (except, in scheme C, for very small or 
very large values of x). Therefore, we conclude that the ambiguity related to subleading 
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terms is substantially reduced if massive terms are treated to (D(a 2 ), i.e. using schemes 
B or C. If massive terms are treated to 0(a s ), i.e. using scheme A, our results suggests 
that the accuracy can be improved using a threshold suppression, and in particular using 
a damping factor to treat threshold effects. It is interesting to observe that, for F2, right 
at threshold the difference between scheme A and scheme B appears to be non-negligible, 
i.e. 0(a 2 ) have a significant impact at threshold, rather larger than the typical size of the 
ambiguity related to threshold terms. 

In the case of Fl results at threshold in scheme A show numerical instabilities due to 
imperfect cancellation between the ZM and MO contributions, which as shown in Fig. [10] 
are much larger than the massive result. They are the consequence of the fact that the 
massless approximation is completely inadequate at threshold for this observable. As soon 
as the scale is raised above threshold this may lead to very large unphysical contributions 
even in the FONLL scheme unless a strong threshold suppression is adopted, such as 
that provided by a damping factor. This is even more the case in scheme C: due to 
the instability displayed in Fig. [10] results without threshold treatment, and even those 
obtained using x scaling are unphysical except at very large x. 

4.3 Comparison of results at different orders 



Finally, in Figs. l20M21l we compare results for the structure functions F2 and Fl determined 
at various perturbative orders, i.e. in the three schemes A,B,C. All results are now deter- 
mined using our default option, namely, the FONLL method with the threshold treated 
using a damping factor: the curves are thus the same as those shown in Figs. 11811 191 for the 
damping factor prescription. For comparison, we also show the simple massless-scheme 
result Eq. (P), determined at NLO (0(a)). 

First, it appears that effects due to the heavy quark mass are not only large in the 
threshold region, but in fact still sizable at Q 2 ~ 10 GeV 2 for F%, and for Fl quite large 
even at Q 2 ~ 100 GeV 2 . Also, it is clear that while the B-C scheme curves are very close 
to each other, the curve in the A scheme differs substantially from them. This means that 
scheme B is an improvement over scheme A, in that it already includes most of the 0(a 2 ) 
effects which are fully included in scheme C. 

Of course, because of asymptotic freedom, differences between results obtained using 
different perturbative schemes become smaller at large Q 2 . Note however that here, for 
illustrative purposes, structure functions have been determined with a fixed set of PDFs. 
In a realistic situation, the structure function would be fixed by the data and the PDFs 
would be fitted. 



5 Mass singularities in F2h 

The heavy contribution to structure functions is experimentally accessible, and indeed 
the experimental determination of the charm and beauty structure functions F% c and 
i<2b has attracted considerable interest recently, in particular at HERA: specifically, the 
kinematic coverage of Fi c data in the (x, Q 2 ) plane for the combined HERA-I dataset [12] is 
shown in Fig. [22l However, the experimental definition of heavy structure functions differs 
somewhat from the definition, given in Sec. 13.21 where Fh was defined as the contribution 
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Figure 16: Same as Fig. [14] but for the 0{a 2 s ) contribution to F2 C (x,Q 2 ). 
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Figure 17: Same as Fig. [16] but at large x. 

to the structure function F obtained when only the electric charge of the heavy quark 
is nonzero. Rather, the experimentally measured heavy quark structure function, which 
we will denote by FX, is defined as the contribution to the structure function F from all 
processes in which there is a heavy quark in the final state [43H50] . 



The difference between F^ and FX is potentially significant, because F^ is affected 
by mass singularities in the limit m — > 0, due to the fact that heavy-quark production 
contributions in which the virtual photon couples to a light quark such as shown in Fig. [2] 
are included in FX, but virtual corrections such as shown in Fig. [5] are not. In contrast, 
using the alternative definition FX neither contribution is included, and both contribute 
to Fi, leading to a cancellation of potential mass singularities. 

For finite values of the heavy quark mass these mass-singular contributions to FX are 
of course finite, but enhanced by double logs (powers of L 2 ): the first diagram of Fig. [2] 
leads to a contribution of order a 2 log 3 Q 2 /m 2 [3"2|l37]. The three logarithmic powers have 
the following origin: one arises from the collinear singularity in the emission of the gluon 



37 




Figure 18: The FONLL expression for F 2c in the three schemes for perturbative ordering 
A-C of Sec. 13.31 and either with no threshold suppression terms, or with the damping 
factor or \ scaling suppression at threshold. The structure functions are plotted as a 
function of x for fixed Q 2 = 2, 4, 10, 100 GeV 2 (from bottom to top). 



from the light quark, one from the collinear singularity of the gluon splitting into the 
hh pair, and one is due to the gluon becoming soft. The latter log arises because the 
contribution to the coefficient function from the diagram of Fig. [2] is singular at z = 1: 
the convolution integral with the PDF up to the kinematic limit z = Q 2 /(Q 2 + m 2 ) < 1 
then leads to an extra log whatever the behaviour of the PDF. At higher perturbative 
orders Fh then receives double-logarithmic contributions of the form a 2+k L 3+2k . These 
contributions could in principle be controlled experimentally by tagging both the heavy 
quark and antiquark and introducing a cutoff on the invariant mass of the heavy quark- 
antiquar pair [37] . 

In conventional parton fits, Fh is usually computed and compared to data, even though 
the data really refer to F^. Furthermore, even if the theoretical expression F^ were imple- 
mented in a parton fit, one may still wonder whether this quantity may be subject to large 
higher-order corrections, because of the aforementioned double logs. It is thus interesting 
to assess the size of the difference between F^ and F^ both at lowest nontrivial order and 
at higher orders. 



5.1 Order a 2 

We have computed the ratio 

M 2c {x,Q) = (83) 

Fi c "(x,Q 2 ) 

for the more practically relevant case of the charm structure function F 2 . We use the 0(a 2 s ) 
expressions for F 2 (n c l \x, Q 2 ) and F 2 {n c l \x, Q 2 ) of Sec. I3T21 the numerator of Eq. flS3]) then 
is the contribution from the diagram of Fig. [2j Since we are mostly interested in the 
behaviour near threshold, the massive scheme is adopted throughout. The conventional 
Les Houches PDFs of Sec. H] are used. Results are shown in Fig. [23l For comparison, we 
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Figure 19: Same as Fig. [18] but for the structure function F^ c . 



also show in Fig. [23] the relative size 
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of the contribution from the light-quark initiated diagrams 



'J i .</; 



(84) 



(shown in Fig. [3]) to 



the heavy structure function f!^. It appears that the light-quark initiated contribution 

to i*2e °^ though moderate, can amount to as much as 10% of Fic m the HERA 

region and it is thus rather larger than the mass-singular contribution Fig. [23l which 
remains at the percent level in this region. Therefore, within this kinematics the mismatch 
Eq. (|83|) between the two definitions of the heavy quark structure function is negligible 
even when the measured light-quark initiated countribution (which includes both of the 
contributions of Figs. [23l24|) is not small. 



Even so, one may wonder whether these results still hold once experimental cuts are 
accounted for. To this purpose, we have computed 

Aa vis , c == ?vis ' c ~ ^ , (85) 

"Vis,c 

where cr v i SiC is the so-called "visibile" reduced cross section (after subtraction of the con- 
tribution from Fl), i.e. the contribution to the DIS charm production cross section with 
experimental cuts. We have used the HVQDIS Monte Carlo program |51j and massive- 
scheme ZEUS-S PDFs |52j : hadronization and charm decay corrections have been ne- 
glected, though their inclusion would be straightforward. In order to define the visible 
cross section, we have assumed > 1.5 GeV and —1.6 < rj c < 2.3, which correspond 
to the acceptance of the recent ZEUS muon analysis [19] when hadronization and charm 
decay corrections are neglected. 

Our results are summarized in Table [lj where we also tabulate the percentage ratio 
Eq. (j83|) already shown in Fig. [23] but now also computed using HVQDIS and ZEUS-S 
PDFs. We have adopted an (x,Q 2 ) binning similar to that of the upcoming combined 
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HERA F% c dataset. We observe that, in qualitative agreement with the results of Fig. [23l 
the percentage discrepancy is always < O (1%), reaching > O (2%) only for the largest Q 2 
bins, where statistical uncertainties are anyway much larger. Clearly, with this choice of 
kinematics, results are only marginally affected by experimental acceptances. We conclude 
that in the HERA kinematic range the discrepancy between F^ and F^ at 0(a 2 ) is at the 
level of the percent, even with experimental cuts. It is interesting to observe that at small 
x and not too high Q 2 this contribution grows, and it should then be properly accounted 
for if F2 c were measured at a future higher energy electron-proton collider [53] . 

5.2 All— order resummation 

When more gluon splitting processes are inserted in the gluon propagator in Fig. [21 before 
the final splitting into the hh pair, they lead to corrections enhanced by further powers 
of a s log 2 Q 2 /m 2 (i.e. double log enhanced). The effect of the complete resummation of 
these double logarithms was studied in detail in Ref. [54] (see also Ref. [55]), in the case 
of heavy flavour production associated with the production of a gluon jet at a scale Q 2 . It 
amounts to an enhancement of the production cross section by a factor n g (Q 2 , K 2 ), where 
K 2 is the virtuality of the hh pair. The way n g is determined and its precise definition 
are discussed Ref. [53]; here, we can interpret it as the enhancement of the light-quark 
initiated heavy quark contribution of Fig. [21 i.e. the numerator of Eq. (j83H . whose relative 
size is shown in Fig. [23l We have K 2 > Am 2 , and Q 2 can be identified with the DIS Q 2 
scale. 

The reason why the scale Q 2 should be chosen as the scale for the evaluation of n g , as 
opposed to W 2 , deserves a comment, given that in the small x limit these two scales are 
widely different. In fact, if we assume a behaviour like l/x 1+<5 for the light quark density 
in the small x limit, the contribution of the first diagram of Fig. [2] has the form 

l f (£) L^ 2 \z,Q 2 /m 2 ) q (x/z) « q(x) (J) jf dzz S L™ + 2 \z, Q 2 /m 2 ) . 

(86) 

Since L NS >( 2 ) is a non-singlet coefficient function, its z integral is finite at the lower limit, 
and the result grows like log 3 Q 2 /m 2 (see Eq. (D.8) and subsequent comment in Ref. |32j . 
see also Ref. [37]). This confirms the scale choice of Q 2 rather than W 2 . 

The function n g is plotted in Fig. [25] as a function of K 2 for different scales Q 2 : at a 
scale K 2 = 4m 2 = 8 GeV 2 , we see that the enhancement is below 20% for Q 2 < 100 GeV , 
and below 90% for Q 2 < 1000 GeV 2 . Hence, the effect is moderate in the HERA region, 
where the contribution which is thus enhanced amounts to a few percent of the structure 
function in the first place, as shown in Fig. [23j However, the effect becomes large as Q 2 
increases and it could be a significant correction at a higher energy electron-proton collider 
such as the LHeC [53]. Also, in the very small x limit, single log enhanced contributions 
of the form a s logW 2 /m 2 will arise |41j . and eventually prevail on the double logarithms 
discussed above. The effect of these small- x logarithms is unlikely to be important the 
HERA energy regime, but it may deserve further studies. 
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Table 1: The percentage contribution of non- vanishing terms in the e c = limit to F2 C , Eq. (|83l) . 
and to cr^ s e °, Eq. (JHSJ. The grid in (x, Q 2 ) assumed is similar to that of the combined HERA Fz c 
dataset, see Fig. [221 
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Figure 20: The FONLL expression for the structure function i^c with the threshold 
treated using the damping factor Eq. (|18p . in the three schemes for perturbative ordering 
A-C of Sec. 13.31 The structure functions are plotted as a function of x for fixed Q 2 = 
2, 4, 10, 100 GeV 2 (from top to bottom). Results in the NLO massless scheme result 
(same as Figs.[8]l9]) are also shown for comparison (dotdashed). 




Figure 21: Same as Fig. [2T1 but for the structure function F^. 
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Kinematical coverage of combined HERA-I F° dataset 
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Figure 22: Kinematical coverage of F 2c in the combined HERA-I [32] dataset. 



6 Conclusion and outlook 

We have presented a study of the theory and the phenomenology of the inclusion of heavy 
quark corrections in deep-inelastic structure functions, using the FONLL approach that 
had been previously proposed in the context of heavy quark photo — and hadroproduction. 
This approach is suitable for the combination of fixed order heavy quark emission terms 
with the all-order resummation of collinear logs which appear at scales much larger than 
the heavy quark mass. A significant feature of the method is that the perturbative order 
at which the fixed-order and resummed results are obtained can be chosen independently 
of each other in the most suitable way: in fact, we have explicitly considered two different 
NLO implementations (denoted as scheme A and scheme B) in which fixed order results 
of order a s or a 2 have been combined with NLO parton distributions. 

After discussing in detail the method and its implementation to 0(pq), and verifying 
explicitly its consistency, we have studied the impact of heavy quark corrections and their 
ambiguity on the F2 and Fl structure functions. We have found that charm mass effects 
have a significant impact, at the level of 10% on the charm structure function F^ c (for fixed 
PDF) for scales as large as Q 2 ~ 10-20 GeV 2 . The effect is rather larger in the threshold 
region, and also for the Fl c structure function, for which it is a sizable correction even at 
Q 2 ps 100 GeV 2 . For scales Q 2 w 4 GeV 2 there is an ambiguity due to subleading terms 
which are not accurate as threshold, which for F<i c at 0(a s ) is almost as large as the whole 
heavy quark correction. 

We have seen that introducing a suppression term near threshold for these subleading 
contributions (such as provided by a damping factor or by so-called x _scann g) somewhat 
reduces this ambiguity. Comparison to exact 0(a 2 ) results suggests that this threshold 
suppression improves the accuracy of the calculation, though beyond these order there are 
no exact results to compare to. The ambiguity is substantially reduced if 0(a 2 ) heavy 
quark production terms are used within the NLO computation. 

We have finally discussed issues related to mass singularities in the experimental defini- 
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quark initiated term - light quark coupled to photon 




x 



Figure 23: The AF^" percentage ratio Eq. (|83|). in the (x, Q 2 ) plane. The region bounded by 
the dashed lines is the HERA kinematic range of Fig. 



tion of the heavy quark structure functions, which differs somewhat from the theoretically 
most natural definition: we have seen that the impact of these corrections is small in the 
HERA kinematic region, but could become relevant at higher energies. 

Besides giving a simple, flexible and practically viable implementation of heavy quark 
effects, our results provide a framework for the understanding of the impact of heavy quark 
effects on the determination of parton distributions and of the dependence of results on 
details of the procedure due to its theoretical uncertainties. 
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A Appendix: Implementation of the FONLL scheme for 

F 2 {x,Q 2 ) up to G(a 2 s ). 

In this appendix, we collect the relevant explicit formulae for the practical computation 
of the the DIS structure function F2(x,Q 2 ) within the FONLL approach up to O (a 2 ), 
which has been already discussed more formally in Sec. 13.2.21 Formulae for Fl(x,Q 2 ) 
are simpler and can be obtained in an analogous way, and thus are not discussed here. 
Unlike the formulae found in Sec. [3j given for the purpose of illustration, here we focus 
on practical implementation issues, i.e., we state which equations and computer codes we 
have used to arrive to the numerical results for the FONLL structure functions that have 
been presented in Sec. [H 

For simplicity, we denote in this Appendix F<i by F and similarly for all the various 
coefficients, and we write throughout 

a s = 4 n < +1 \Q 2 ), h{y,Q 2 ) = f\ ni+l \y,Q 2 ) . (87) 



As has been discussed in Sec. 13.21 in order to obtain the full structure function F we 
compute separately the "light" F\ and "heavy" F^ contributions. For F^, up to order 
(D(a 2 ), the relevant equations are 



if ONLL (*,Q 2 ) = if } (*,Q 2 ) + it'W), 

if>(x,Q 2 ) = F^ +1 \x,Q 2 )-F^°\ X ,Q% (89) 

i=g,q,q Jx x U KU J 

F^\x,Q 2 ) = f^C^U^L^My^) (91) 



i=9,Q,g 



F^\x,Q 2 ) = x £ f^fe^M 2 ), (92) 

_ , -rJx y \y j 

i=g,q,q,h,n 



where x x = x{\ + 4m 2 /Q 2 ) and are the massless limit of the massive coefficient 

functions Note that in Eqs. (|88H92p the PDFs and the strong coupling are always 

given in the n\ + 1 flavour scheme. 

Up to 0(a 2 ), the massive coefficient functions C^' are given by 



m 2 J (2vr) 2 m 2 



M ( z Q 



m- 

2 V Q 2 
i > \ m z I rrv 



(93) 



The massive coefficient functions c^ g , c^ g , <4 q , c£ q were first computed in [21], together 
with the analog expressions for F^. They are given as phase-space integrals of differential 
cross sections; the corresponding op , cp coefficients (as defined in Ref. [21]) are plotted 
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quark initiated term - heavy quark coupled to photon 




x 

Figure 24: Same as Fig. [23l but for the percentage contribution i?u g ht of light-quark initiated 
terms to F 2c Eq. (JSJ). 




Figure 25: The resummation factor n g (Q 2 , K 2 ) as a function of K 2 , for two values of Q 2 . 
Upper lines have A = 300 MeV, lower lines have A = 200 MeV. 
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in Figs. 6-10 of that reference (the quark coefficient has been subsequently corrected as 
Fig. 6 of Ref. [51]). For the FONLL phenomenological studies of Sec. 0] we have used 
direct numerical integration of the expressions of Ref. [21] , benchmarked against the plots 
in Refs. [21] and [56]; we have not relied upon the interpolation provided in Ref. [57] . 
The massless-limit expressions of the C^' 1 coefficients, needed to compute C^ ! '° were 
taken from a Fortran implementation [5B] of the results of Ref. [32]; they were tested at 
large values of Q 2 against our massive expressions. The MS cross sections i^ ni+1 ^ at the 
NNLO order were computed using the QCDNUM program |39j . 

For the "light" contribution F\ instead, up to order 0(a 2 ), it is convenient to rear- 
range the expression Eq. (|65|) of Sec. 13.2.21 exploiting the fact that all contributions to 
p( d ) Eq_ (|67p cancel except the heavy quark contribution, and that the massive-scheme 
contribution with ni flavours is most easily obtained by writing it in terms of the corre- 
sponding contribution with n/ + 1, which can then be obtained from the QCDNUM code. We 
thus wind up with the following expression: 



F 



FONLL 



(x,Q 2 



F?\x,Q 2 



(x,Q 2 



F 



F 



F 



(d) , 
T 

£ 

i=h,h 

K+i) 



1 dVgtrn+l) ( X 



y 



i-l 



+ x 



ij^h,h,g 



x J2 

i=h,h 

l dy 



D; 



y 

1 dVMni+1) 

x Q 2 



y 



y m z 



a. 



My,Q 2 



(94) 
(95) 

(96) 



where we have defined 



(£)' 



4 2) *A )-s(i-z) 



dzL q 2) 



(z) 



dni 



(nj+l) 



(97) 



(2) 

where the function L q is given in formula (A.l) of Ref. [32], with the further assumption 

Q 2 \ _ . 1 



L (2) 
1 



m 



2 ) O^s 







for z > 



1 + 



4m 2 



(98) 



The function K qq is given by 



K qq (z,L) 
K qq {z) 



K qq (z) + 
'14 



L 2 2T R 



CfTr 
CfTr 



2 3 ~ qq 

2 n 



P^°(z)-LA qq (z), 



1 - z 
l + z 1 



28 



-log z + -\ogz + - 



+ (1-*) 



2 1 W \ 4 / n 

3 log, + y +-(!-,) 



2 , 13 

— log z -\ 

3 S 9 



(99) 

(190) 

+ 

(101) 
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and the NLO quark coefficient function (see for example Ref. [59]) is 



l + z 2 ( 1 -z 3\ I, 



(102) 



The term proportional to 5(1 — z) in Eq. (j97|) is the heavy flavour virtual correction, and 
it is needed to enforce the Adler sum rule, which thus determines its coefficient. We can 
easily see that 

Di,i (-,^,a s ) = B U { X -^ a \-cf^ (~,a s ) , (103) 
\y m 2 J \y m 2 J \y J 

so that, using also Eq. (f6"3"jh we have 

F? l) = x £ f^By (-X 2 ,a s ) fi(y,Q 2 )- (104) 

.rfrJx y \y m 2 j 

In Eq. (jnSJ) , the term F^ ni+l ^ is evaluated directly using the QCDNUM program. All 
remaining terms must be evaluated independently. Note that the second term on the 
r.h.s of Eq. (|96|) is identical to Eq. (J95j) : these terms thus would cancel exactly, but 
the cancellation does not happen if one uses a threshold prescription such as Eq. (|19p or 
Eq. (|20p . because the threshold prescription only affects but not . The coefficients 
of ni in C^ H+1 ^ and C^J^ have been evaluated using the same equations as the QCDNUM 
package, given in Refs. [60] and [61], since they have to compensate the QCDNUM result in 
p(ni+i) ^ N ^ ce that, according to Eq. (|7ip Dn vanishes in the massless limit. 

All the 0(a 2 s ) results given so far, used with NNLO PDFs, correspond to the scheme 
C of Sec. 13.31 The 0(a s ) results (scheme A) were obtained using the same formulae, trun- 
cated at order a s , and using NLO evolution and massless coefficient functions in QCDNUM. 
The alternative NLO implementation, denoted by scheme B in Sec. I3.3[ uses the NLO 
approximation in the massless result, and the full 0(a 2 s ) result in the massive one is also 
obtained in a similar way. The only caveat is that in the subtracted term in Eq. (|9ip we 
replace the massless limit of the C^' 1 (denoted by C^' ^' 1 ) with 

<% m (*, S) - cJ3** {*%)- <% m <M). (105) 

i.e. the non-logarithmic term is excluded from the subtraction, according to the discussion 
in Sec. I 
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